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Working Paper: Exact Convex Relaxation for 
Optimal Power Flow in Tree Networks 

Lingwen Gan, Ufuk Topcu, Na Li, Steven Low 
Abstract 

The optimal power flow problem is non-convex. We modify its voltage constraint to be "slightly" more conser- 
vative, and then relax the modified problem to a convex problem in tree networks. A sufficient condition for exact 
relaxation is provided. This condition can be checked before solving the relaxation, and holds in various test networks, 
including IEEE 13, 34, 37, 123-bus test networks and a real distribution network with high penetration of distributed 
generation. 

I. Introduction 

The optimal power flow (OPF) problem seeks to minimize a certain objective function, such as power loss and 
generation cost, subject to physical constraints including Kirchoff s laws, voltage regulation constraints, and power 
injection constraints. There has been extensive research on OPF since Carpentier's first formulation in 1962 fT\, 
and surveys can be found in ||2l, lEl. EIj 0, O- OPF is in general non-convex, and a lot of algorithms have 
been proposed to approximate OPF or relax OPF to convex problems. Relaxation methods have the potential of 
providing exact OPF solutions, and are the focus of this work. We call a relaxation exact if every solution to the 
relaxed problem also solves the original problem. 

In transmission networks, which are usually mesh networks, semi-definite relaxation (SDR) has been proposed to 
solve OPF 121, m, ini, but whether or when SDR is exact can only be checked after solving SDR. In distribution 
networks, which are usually tree networks, different convex relaxations for OPF have been proposed based on the 
DistFlow power flow model proposed in ifTOl . ifTTI . Reference lfT2l puts forward a second-order-cone relaxation 
(SOCR) for OPF, and proves that SOCR is exact if there are no lower bounds on power injection. This condition 
. can be checked before solving SOCR. See also lfT3l . lfT4l for other relaxations for OPF. 
^^^^^^ The lower bounds on power injection are important for various applications including demand response ifTSll 
and VoltA^AR control lfT6l . Motivated by this shortcoming of (12], references ifTTl . ifTSl work on other sufficient 
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conditions for the exactness of SOCR. Though references ifTTll . ifTSl bring back the lower bounds on power injection, 
they remove the upper bounds on vohage. The upper bounds on vohage are important for vohage regulation, 
especially when there are distributed energy sources increasing the likelihood of over-voltage, i.e., voltage exceeds 
voltage upper bound. 

Contributions of this work include the following. Firstly, we show that SOCR is not always exact, especially 
in the presence of distributed generation. Furthermore, over-voltage is a cause for inexact SOCR. Secondly, we 
modify OPF by using a "slightly" more conservative voltage constraint to deal with the over-voltage issue. The 
modified OPF (MOPF) problem is still non-convex, therefore we relax it to a convex problem called relaxed MOPF 
(RMOPF). Thirdly, we provide a sufficient condition for the exactness of RMOPF. This condition can be checked 
before solving RMOPF, and holds in various test networks including IEEE 13, 34, 37, 123-bus test networks 
|fT9l and a practical network of Southern California Edison ll20l with high penetration of distributed generation. 
The condition roughly imposes that reverse power flow — power flowing from the loads to the substation — in the 
network should not be too large. At last, we present our proof technique. While prior works ifTTll . ifTSi focus on the 
dual method — prove properties of the dual variables — to derive sufficient conditions for exact relaxation, this work 
explores the primal problem directly. We show that the results in IfTTll can be re-established using the technique in 
this paper. f 

The rest of the paper is organized as follows. In Section |ll] we introduce the OPF problem and its relaxation 



SOCR, as well as a counter example where SOCR is not exact. In Section III we propose the modified OPF 
problem MOPF, give a convex relaxation RMOPF for MOPF, provide a sufficient condition for the exactness of 
RMOPF, and check that this condition holds in various test networks. In Section IV we present the proof technique, 
and use this technique to re-establish the results in I fTl. 



II. The Optimal Power Flow Problem 

A. The Optimal Power Flow Problem jr 

Consider a tree distribution networl|^ that consists of iV + 1 buses. Index the substation bus by and the other 
branch buses by i = 1, . . . , A^. Let A/^ := {0, . . . , A^} denote the set of all buses. If two buses i,i ^ M are connected 
by a line, we denote i ^ j; otherwise, we denote i oo j. 
^ Definition 1: The ordered set {io,ii, ■ ■ ■ ,ik) ^ Af where fc G N is called a path from bus io to bus ik, if 

it ~ it+i fort = 0,...,k-l. 

There exists a unique path between two arbitrary buses in a tree ||2T1 . and we denote the unique path from bus 
to bus i by Vi, for i e J\f. Each transmission line in the network connects two buses, say i and j. Then, either 
i £ Vj or j E Vi, i.e., either bus i is on the path from bus to bus j, or vice versa. We denote the transmission 
line by if i E Vj, and {j, i) otherwise, i.e., we put the bus that is closer to bus in front. Let 

^ {{hj) I i ~ j and i e Vj} 



Most distribution networks ai'e tree networks. 
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denote the set of all transmission lines. 

For each bus i G J\f, let Pi and qi denote its real and reactive power injection respectively, let Vi denote its 
voltage, and define Vi := \ Vi\^. Bus — the substation — ^has fixed voltage magnitude, i.e., vo is a constant. For each 
transmission line G C, let Vij and Xij denote its resistance and reactance respectively, let Pji and Qji denote 
the real and reactive power flow from j to i respectively, let 7^ denote the current flowing from i to j, andi 



:= |7jj|^. Define 



X := {Xij,{i,j) e Cf 
P:= {Pji,{i,j)&Cf 



q := (gi,...,gjv)^ 
V := {vi, . . .,Vn)'^ 



idjdRme 

2> 



i.e., use the letter without subscript to denote a column vector of the corresponding quantity. For any two vectors 
a = (ai, . . . , Qk) e K'^ and b ■ 



. . . , 6fc) e M'' where fc e {1, 2, . . .}, define the relations ^, ^ as 

ay b Oj > 6j for i = 1, 
a^b 4^ ai> b j fg^i = 1, 
a ^b Oj <mi for^|= 1, . . . , fc, 

a^b 4^ Oi < 6i for i = 1, . fc. 



4^ 



i - j 
i oo j 

r^ 



p 

Qji 
lij 
f- ■ 



the set of buse 
bus i andgiffl-e cormected 



b'l^ and'j are not connected 
the path from bus to bus i 
a transmission line, i ^ j and i £ Vj 
the set of transmission lines 



an^iP^^i 



real power injection at bus i 
reactive power injection at bus i 
voltage at bus i 



resistance of (z, j) 
reactance of 

real power flow from j to i 
reactive power flow from j to i 
current flowing from i to j 

I - = 



The optimal power flow (OFF) problem is 
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OFF: 



s.t. 



p,q,P,Q,e,v 

Vi — Vj 'Zl^Tij Pj'i -\- 

where r, x, Vj, u, p, p, q, q, and vq are given constants. 

Remark 1: The objective function in ([T]) is the power loss, and can be generalized to any function of the type 




/ '^^i^'-^J +9{p,q), 

where / : M — > M is strictly increasing and g : x — > M is arbitrary, without changing the results in this 
paper. For ease of presentation, we work with the objective function in ([TJ. 

Remark 2: Equations (|2|-(|5]) are due to the underlying physical laws that power flow satisfies lITOll . ifTTl . and 
called the DistFlow model in the literature. It is claimed in ll22l that once substation voltage wq and power injection 
p and q are specified, then there exists a unique practical P, Q, v satisfying the DistFlow model (|2])-(|5]). The 
control is to find the optimal power injection p and q, and let the other variables P, Q, t and v be determined by 
the DistFlow model (l2])-(|5]l. ^ 

Remark 3: Equation ([6| is the voltage regulation constraint. Voltage is usually expressed in terms of the per 
unit value. For instance, the voltage constraint can be 0.95p.u. < \yi\ < 1.05p.u. |23|, which translates into 

Remark 4: Equation (|7| are the power injection constraints. In applications including demand response and 
VoltA'AR control, power injection constraints are usually considered to be box constraints as (|7]i, but results in this 
paper extend to arbitrary power injection constraints (they do not even need to be convex). 

B. The Second-Order-Cone Relaxation 

OPF is non-convex due to the non-affine equality constraint in (|5]l. An approach (see lfT2l ') is relaxing this 
constraint to the inequality constraint 

P2. + 02. 

% > ^' , {hj)&C. (8) 



August 21, 2012 



DRAFT 



5 



Equation ([8]) is equivalent to a second-order-cone constraint. Therefore, we call the following convex problem the 
second-order-cone relaxation (SOCR) of OPF. 



SOCR: 



S.t. 



E 

p,q,P,Q,£,v 



If (|5]l is satisfied at an optimal solution w°p^ to SOCR, then w°p*' also solves OPF. 

Definition 2: A relaxation is exact if every of its optimal solutions solves the original problem. 
When a relaxation is exact, we can solve the original problem by solving the relaxed problem. The original problem 
is usually non-convex and difficult to solve, while the relaxed problem is usually convex and easy to solve. 

In this section, OPF is the non-convex original problem, and SOCR is the convex relaxed problem. SOCR is 
exact if and only if every of its optimal solutions satisfies the constraint in (|5]l. It is proved in lfT2l that if there 
are no lower bounds p and q in (|7]i, then SOCR is exact. In applications including demand response and VoltA'AR 
control, lower bounds p and q exist, requiring further analysis. 



C. SOCR is Not Always Exact. * 

SOCR is not always exact. Furthermore, when it is not exact, the solution to SOCR is physically meaningless. 
We illustrate this through a 2-bus network in Fig. [T] Bus is the substation and has fixed voltage magnitude = 1- 



V 



Fig. 1. A pathological two-bus circuit. 



The branch bus has renewables generating Ip.u. real power, of which p E [0, 1] will be injected to the grid and the 
rest I — p will be curtailed. The reactive power injection q is 0. Line resistance and reactance are r = a; = 0.1, 
and the voltage constraint is 0.9 < < l-l- The objective is to minimize the sum of power loss r£ — O.ll and 
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curtailment 1 — p. The OPF can be written as 

min 0.1£+{l~p) 
over p, q, £, v 



s.t. 



V = 1 + 0.2p + 0.2g- 0.02£, 

p2 ^ q2 



and the corresponding SOCR is 



over 



i = 

V 

0.9<v< 1.1, 
< P < 1, q = 

0.1£+{l-p) 
p,qj,v 



S.t. u = 1 + 0.2p + 0.2g- 0.02£, 




e > 



p + q 



0.9<v < 1.1, 



The solution to ([T0| i^ 



0<p< 1^ = 0. 
0.7705, g = 0, ^ = 2.7049, v = 1.1, 



which violates £ — {p^ + q^)/v, therefore it is physically meaningless. 

There is a tendency to operate at p = 1 since curtailment is penalized, but large p leads to over-voltage, i.e., 
V > 1.1. The real optimal solution for (|9]l — if the power loss term 0.1^ is negligible — is to increase p up to p°P* 
where v = 1.1. However, after relaxing (|9| to ( [T0] i, it is possible to increase p beyond p°P' to obtain a smaller 
objective, while allowing £ > (p^ + q'^)/v to satisfy v < 1.1. 

To summarize, SOCR is not always exact, and attains a physically meaningless solution when it is not exact. 
Furthermore, over-voltage is a cause of inexact relaxation. 



III. A Modified OPF Problem 
In this section, we modify OPF to deal with the over-voltage issue, by changing the constraint u < iJ in (|6| to 



,,lin 



< V, where v is a linear approximation as well as an upper bound for v. We first propose a modified OPF 



problem MOPF in Section III-A and then compare MOPF with OPF in Section III-B and III-C Afterwards we 



provide a convex relaxation RMOPF for MOPF as well as a sufficient condition for the exactness of RMOPF in 



Section |III-D| This condition can be checked before solving RMOPF, and holds in various test networks as shown 
in Section HiTeI 



-We solve SOCR-example using the Matlab toolbox CVX [24l. 
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A. A Modified OPF Problem 

We modify OPF to deal with the over-voltage issue, which may cause inexact relaxation (see the example in 



Section II-C i. The way we modify OPF is to replace the constraint w < t; in (|6]l by a new constraint w''" < v. Here 
w''" is a linear approximation as well as an upper bound for v. ^m.. 

We start with introducing the linear approximations P''", Q''" and to P, Q and v. They are used to study 
optimal placement and sizing of shunt capacitors in ITOl . ifTTl . to minimize power loss and balance load in 125], 
and to control reactive power injection for voltage regulation in lfT6l . The linear approximations P'' 
are solution to 

(j,fc)e£ 

Qji ^ 1j ^ Qkji 
Vi — Vj liVijPjj^ ~\- XijQji)^ 




e., ignoring the I terms in (|2]l-(|4|. The ^ terms are much smaller than the other terms for a broad class of distribution 



networks including the ones considered in this paper. Linear approximations P'"\ Q''" are valid when rijtij <C |Pj 



lin ie 



and Xijlij <C \Qji\ for G C, i.e., power loss is much smaller than power flow. Linear approximation v 
valid when we further have r^j \Pji{p^q)\ <^ Vj and Xij \Qji{p,q)\ ^ Vj for G C, which happens to be the 
requirement of voltage stability lfT6l . Therefore, w''" is valid when power loss is insignificant and voltage is stable. 
The approximations have closed-form expressions '"J'^ 

Pft ^ E (14) 
E ^fe, («,j)e A (15) 




(16) 

Note that given p, q, vq, the quantities P''", Q''" and u''" can be computed without solving the nonlinear power 
flow equations ([2|-(|5]l. 

Equations ([T4]i and ( [T5] l imply that Pjj" and Q^'" are equal to the sum of their subsequent real and reactive power 
injections respectively for G C. It follows that P''" is a linear function of p, Q''" is a linear function of q, 
and w^'" is a linear function of p, q. Furthermore, the following lemma shows that P''", Q''", and u''" are upper 
bounds on P, Q, and v. 

Lemma 1: Let {p,q,P,Q,i,v) satisfy (|2|-(|4| and (|8]l, then 

P ^ P""(p), Q ^ Q'jfiq), and v ^ ^""(p,?). 
Proof: In Appendix [A] ■ 
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We replace the voltage constraint v < v in (|6]l on i; by the voltage constraint w''" {p, q) < v in ( [T7] l on 
to obtain the modified OPF (MOPF) problem 
MOPF: 



min 








over 


p,q,P,Q,£,v 


s.t. 


0; 




v^v, (?) ^ V. 



^^^^ 



B. Comparison of MOPF and OPF 

Problem MOPF and OPF only differ in the voltage constraint, i.e., the upper and lower bounds on the voltage 
profile. Consider the hypothetical scenario in Fig.|2] The dashed curves illustrate the voltage upper and lower bounds 
in OPF. MOPF has the same voltage lower bound as OPF. Since v < v^™, the constraint t;''" ^ u in MOPF can be 
interpreted as having a stricter upper bound on the voltage, as illustrated by the dotted curve. Since MOPF has a 
stricter voltage constraint than OPF, every feasible voltage profile for MOPF is also feasible for OPF, but voltage 
profiles that are feasible for OPF are not necessarily feasible for MOPF. For instance, in Fig. |2] voltage profile B, 
shown by a dash dot curve, is feasible for both MOPF and OPF, while voltage profile A, shown by a solid curve, is 
feasible only for OPF. On the other hand, when power loss is insignificant and voltage is stable — normal operation 
points satisfy these, d''" and v are "close," therefore the two upper bounds in Fig. |2] are "close." 



Voltage 



> (P-u.) 




1.05 



0.95 




Bus 



Path -P, 



Bus j 



Fig. 2. Voltage constraints for MOPF and OPF along a path Vj for some j S A^. Dashed curves illustrate the voltage upper and lower bounds 
in OPF. MOPF has the same voltage lower bound, but a smaller voltage upper bound illustrated by the dotted curve. Both voltage profiles A 
and B satisfy the voltage constraint in OPF. but only B satisfies the voltage constraint in MOPF. 



We now formalize the intuition from Fig. |2] Let J^opf denote the feasible set for OPF, and Jiviopf denote the 
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feasible set of MOPF, then Ji^opf ^ -^opf- Define 

e := vaax v^™{p,q) ~ Vi{p,q) 
s.t. P ^ P ^P, 

1^1,. ..,N 

as the maximum deviation from q) to v{p,q), then 

Mp, q)<Vi-e)^ q) < Vi) 

Consider OPF with a stricter voltage upper bound Vi < Vi — e 
OPF-e: 

min ^ njiij 

over p,q,P,Q,£,v 
s.t. (|2]l-(l5]l, Q; 

< Vi <Vi - e, i = 1,. . . ,N. 




It follows from ^TE) that J"opf-£ ^ JSiopf- Hence, 

^ -7^0PF-e C J"mOPF Q J'OPF, 

i.e., MOPF is "sandwiched" between OPF and OPF with stricter voltage upper bound. 




C. Evaluate e 

As seen in OPF-e, the maximum deviation e from ti''"(p, q) to v{p, q) can be considered as the conservatively in 
voltage upper bound. If = 0.9, Vi = 1.1 for alii £ J\f and e = 0.01, then the voltage constraint is 0.9 < Vi < 1.1 
for OPF and 0.9 < < 1.09 for OPF-e. For voltage profiles that satisfy the voltage constraint of OPF but violate 
that of OPF-e, we must have 1.09 < Vi < 1.1 for some i e Af. The voltage Vi is close to though still below the 
voltage upper bound in OPF, and such operating points are to be avoided due to safety consideration. The smaller 
e is, the "closer" is OPF-e to OPF, and the less favorable are the operating points that are feasible for OPF but 
infeasible for OPF-e. 

It has long been accepted in the literature ifTOll . ifTTI . ||25]| . lfT6l that the gap e between and v is negligible 
for practical networks, and we provide an empirical evidence that e = 0.0036 for the network in Fig. |3] It is 
a practical distribution network of Southern California Edison (SCE) |'20'| with high penetration of distributed 
generation (11.3MW peak load and 6.4MW nameplate distributed generation capacity). Line impedances, peak spot 
load, and nameplate ratings of shunt capacitors and photo voltaic generators of the network are summarized in 
Table |l] In evaluating e, we assume that the loads are drawing peak spot apparent power at power factor 0.97, all 
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Fig. 3. An industrial feeder of SCE with high distributed generation (photovoltaics) penetration. Bus 1 is the substation and the 6 loads attached 
to it model the other feeders on this substation. 

^ TABLE I ■ 

Line Impedances, Peak Spot Load, and Nameplate Ratings of Capacitors and PV Generations of the Network Shown in 



Network Data 



Line Data 


Line Data 


Line Data 


Load Data 


Load Data 


PV Generators 


From 


To 


R 


X 


From 


To 


R 


X 


From 


To 


R 


X 


Bus 


Peak 


Bus 


Peak 


Bus 


Nameplate 


Bus. 


Bus. 


(Q) 


(Q.) 


Bus. 


Bus. 


m 


(n) 


Bus. 


Bus. 


(Q.) 


(O) 


No. 


MVA 


No. 


MVA 


No. 


Capacity 




























1 


2 


0.259 


0.808 


8 


41 


0.107 


0.031 


21 


22 


0.198 


0.046 


1 


30 


34 


0.2 






2 


13 








8 


35 


0.076 


0.015 


22 


23 








11 


0.67 


36 


0.27 


13 


1.5MW 


2 


3 


0.031 


0.092 


8 


9 


0.031 


0.031 


27 


31 


0.046 


0.015 


12 


0.45 


38 


0.45 


17 


0.4MW 


3 


4 


0.046 


0.092 


9 


10 


0.015 


0.015 


27 


28 


0.107 


0.031 


14 


0.89 


39 


1.34 


19 


1.5 MW 


3 


14 


0.092 


0.031 


9 


42 


0.153 


0.046 


28 


29 


0.107 


0.031 


16 


0.07 


40 


0.13 


23 


1 MW 


3 


15 


0.214 


0.046 


10 


11 


0.107 


0.076 


29 


30 


0.061 


0.015 


18 


0.67 


41 


0.67 


24 


2 MW 


4 


20 


0.336 


0.061 


10 


46 


0.229 


0.122 


32 


33 


0.046 


0.015 


21 


0.45 


42 


0.13 






4 


5 


0.107 


0.183 


11 


47 


0.031 


0.015 


33 


34 


0.031 





22 


2.23 


44 


0.45 


Shunt Capacitors 


5 


26 


0.061 


0.015 


11 


12 


0.076 


0.046 


35 


36 


0.076 


0.015 


25 


0.45 


45 


0.2 


Bus 


Nameplate 


5 


6 


0.015 


0.031 


15 


18 


0.046 


0.015 


35 


37 


0.076 


0.046 


26 


0.2 


46 


0.45 


No. 


Capacity 


6 


27 


0.168 


0.061 


15 


16 


0.107 


0.015 


35 


38 


0.107 


0.015 


28 


0.13 










6 


7 


0.031 


0.046 


16 


17 








42 


43 


0.061 


0.015 


29 


0.13 


Base Voltage (kV) = 12.35 


1 


6000 kVAR 


7 


32 


0.076 


0.015 


18 


19 








43 


44 


0.061 


0.015 


30 


0.2 


Base kVA = 1000 


3 


1200 kVAR 


7 


8 


0.015 


0.015 


20 


21 


0.122 


0.092 


43 


45 


0.061 


0.015 


31 


0.07 


Substation Voltage = 12.35 


37 


1800 kVAR 


8 


40 


0.046 


0.015 


20 


25 


0.214 


0.046 










32 


0.13 






47 


1800 kVAR 


8 


39 


0.244 


0.046 


21 


24 
















33 


0.27 
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shunt capacitors are switched on, and distributed generators are generating real power at their nameplate capacities 
with reactive power (p equals p equals the real power injection in this case, and q equals q equals the reactive 
power injection in this case). We also assume vq — 1. 

The maximum deviation e from w''" to v is considerable only when the power loss is significant with respect to 
the power flow, or the voltage is unstable. Both cases are to be avoided in practice. 



D. Exact Relaxation for MOPF (^^j^ 
MOPF is still non-convex due to the non-affine equality constraint in (|5]), and we relax it to the inequality 

constraint in (|8]l to obtain a convex relaxation RMOPF. 
RMOPF: 



r -f- ■ 



S.t. 



E 

(i,i)e£ 

p,q,P,Q,e,v 



In this section, MOPF is the non-convex original problem, and RMOPF is the convex relaxed problem. RMOPF 
is exact if and only if every of its optimal solutions satisfies the constraint in (|5j. When RMOPF is exact, we can 
solve MOPF by solving RMOPF. ^ ^ 

The main result of this paper is a sufficient condition for the exactness of RMOPF, which can be checked before 
solving RMOPF. Recall the closed-form expressions for P''" and Q''" in ([T4]l-([T5]l, and define 



prev 

ji 



(19) 



for G C as the nonnegative part of Pj"^ and Q^'f- We refer to P"^"^ and Q'^°^ as reverse power flow hereafter. 



Define 



O : 



n 

{k,i)eVi 

E 
E 

{k.,i)eVi 

n 

{k,i)eVi 



1 - 



1 



for i G A^, and note that a^„j for m,n e {1,2} and i £ J\f can be computed before solving RMOPF. The sufficient 
condition for the exactness of RMOPF is 
Theorem 1: RMOPF is exact if 



(20) 



^21 



^22 
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Reader needs to read Section IV-A and IV-B before reading the proof of Theorem [T] in Appendix [P] 

Note that condition pO] ) can be checked before solving RMOPF. It can be written in a simpler (and easier to 
comprehend) form 

^' efe,^.)-^^,^^ V(z,,)e/: (21) 




when > and > for i e Af. In practice, > and > hold for i e Af. Hence, one can check 
( |2T| l instead of ( pO] ) (but equivalently) to guarantee exact RMOPF. 

Note that the interval in pT| ) depends on the reverse power flow. When there is zero reverse power flow, 

b^ — and hi — 00 for i £ J\f. The larger the reverse power flow is, the smaller the intervals (b^,bi) are. Define 
the minimum interval 

as the intersection of all intervals in (|2T}. If the minimum interval {b,b) covers the range of rij/xij, then ( |2T| ) 
holds, which implies exact RMOPF. According to ll26l . the r/x range for practical transmission lines is [0.1, 10], 
and we will refer to this range as the practical range hereafter. 



r the intervals {b^,hi^ 



E. Case Studies ^ 

In this section, we show that the minimum interval 5) covers the practical r/x range [0.1, 10] in various test 
networks, including IEEE 13, 34, 37, 123-bus distribution test networks II19II and a practical distribution network 
of SCE EOl with high penetration of distributed generation, ^jr- ^ 

We calculate the minimum interval (fe, fe) for each of the IEEE test networks, and summarize the results in Table 
[n| It can be seen that the practical range [0.1, 10] is covered by each of the minimal intervals. In our calculations, 
we consider the worst case where there is no load, while all shunt capacitors are switched on. This worst case 
maximizes reverse power flow, and minimizes (5, S). Since the IEEE test networks are unbalanced three-phase 
networks, we assume different phases are decoupled, calculate the minimum interval for each phase, and take their 
intersection as the final minimum interval (6, 5) . 

TABLE II 

Minimum Intervals of Different Test Circuits 




test circuit 




IEEE 13-bus 


(0.0119,00) 


IEEE 34-bus 


(0.0049,00) 


IEEE 37-bus 


(0,00) 


IEEE 123-bus 


(0.0067, 00) 


SCE 47-bus 


(0.0366,15.5739) 
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There is no distributed generation in each of the four IEEE test networks. Therefore the reverse real power flow 
picv _ ^jj^ jj. follows that hi — oo for i E Af. This is why b ~ oo for all four IEEE test networks. Distributed 
generation potentially introduces reverse real power flow, and shrinks the minimum interval (6, &) . 

To explore the effect of distributed generation, we revisit the SCE distribution network in Fig. [3] which has five 
photo voltaic (PV) generators and four shunt capacitor banks generating reverse power flow. Its minimum interval 
is calculated to be (0.0366, 15.5739), which covers the practical range [0.1, 10]. In our calculation, we also consider 
the worst case where there is no load, while shunt capacitors are all switched on and PVs are all generating real 
power at their nameplate capacities. This worst case maximizes reverse power flow, and minimizes (6,6). ^ 

We are conservative in checking ( |2T] l. First of all, instead of checking that (6j,6i) covers rij/xij for each 
(?, j) e C, we check that their intersection (6. 6) covers the practical r/x range [0.1, 10] for all practical transmission 
lines. Secondly, we consider the worst case (which is unlikely to happen) where reverse power flow is maximized 
and therefore (6, 6) is minimized. While (6, 6) covers the practical range [0.1, 10] after these two conservative steps, 
condition pT| ) is more widely satisfied. For instance, if we use the load data in Table |l] to calculate the minimum 
interval of the SCE network in Fig. [s] we obtain (6, 6) = (0.0245, 1209.7), which is much larger than the interval 
(0.0366, 15.5739) obtained assuming the worst case. '^^fcli^ 

To summarize, condition ( |2T| is satisfied in all test networks. Therefore RMOPF is exact in all test networks. 

1 ^ 

IV. Proof of the Main Result 

We have given a sufficient condition for the exactness of RMOPF in Theorem [T] and shown that this condition 
holds in various test networks in Section III-E In this section, we present the proof technique for Theorem [T| as 



well as some other theoretical results. We first present the proof technique for Theorem [T] in Section IV-A| and 



IV-B and then use the technique to re-establish the results in prior work ifTTl in Section |IV-C| At last, we study 



the uniqueness of MOPF solutions in Section |IV-D| 

Prior works BlTI , IITSi use the duality theory [21 , Chapter 5] to derive sufficient conditions for the exactness of 
SOCR. The idea is to prove that the optimal dual variables corresponding to (|8]l are always positive. It follows from 
complementary slackness that the equality in ([8| is attained at any optimal solution to SOCR, therefore SOCR is 
exact. We call this method the dual method. 

^ We explore a different method to derive sufficient conditions for the exactness of RMOPF in this work. Instead 
of working with the dual variables, we directly work with primal variables p, q, P, Q, £, and v. The idea is to 
prove that for any feasible point w = {p,q,P,Q,£,v) of RMOPF that violates (|5]), we can find another feasible 
point w' = {p' , q' , P' ,Q' ,£' ,v') of RMOPF that has a smaller objective value. It follows that the optimal solution 
w°P' to RMOPF must satisfy (|5]), therefore RMOPF is exact. We call this method the primal method. 

A. The Primal Construction 

The key to primal method lies in constructing a new feasible point w' of RMOPF that has a smaller objective 
value, given any feasible point w of RMOF that violates (|5|. For ease of presentation, we present the construction 
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in a one-line network shown in Fig. |4] and the construction in a general tree network can be found in Appendix 
[D| In the one-line network, we abbreviate rij, Xij, Pji, Qji, and iij by rj, Xj, Pj, Qj and £j respectively. With 



I ri,xi 



Vn 



Pq, Qg 



'1 



rN,XN 



Pn, Q 



N, (-N 



Pn^Qn 



Fig. 4. Simplified notation for a one-line network. 

this simplified notation, we re-state MOPF and RMOPF as follows. 
MOPF-line: 



N 



mm ^ rj^i 
1=1 

over p,q,P,Q,i,v 
s.t. Pn =Pn, Qn = Qn 



Pi-i ^ Pi~ nii + pi^i, i = l,...,N; 
Qi-i ^Qi- x,e, + i = l,...,N; 

Vi-i = Vi~ 2(riPi + XiQi) 



, . . . , J. » , 



N: 




, . . . , J. T , 



iV; 



N 

mm ^ rill 
1=1 

over p,q, P,Q,£,v 



s.t. 



p2 + r)2 



2^ 



(22) 
(23) 
(24) 

(25) 
(26) 
(27) 
(28) 
(29) 



(30) 



Given any feasible point w ~ {p, q, P, Q, £, v) of RMOPF-line that violates p6l ), let 

m := min < i > 1 > — 
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denote the first bus that violates p6| ), then m E JV. Pick e > satisfying e < '"^ ~im (such e exists according 
to the definition of m), and construct w'{e) := (p' , q' , P' , Q' , v') according to Algorithm[r| First of all. Algorithm 
[T] keeps p and q unchanged, so that w' still satisfies (|28]l-(|29|. The main construction is for P' , Q' and I', after 
which v' is simply constructed to satisfy ( |25] l. 
In the construction of P' , Q' and I' , we set £' as 



max{Pf,P,^}+max{Q',^Q^} 



& 



Ik k > m, 

^'k = \ ik — ^ k — m, 

k < m. 

That is, we do not change £k for k > m; we reduce £m to £'„^ — £„i — e; and we modify £k for k < m so that 
constraint ( [SO] ) remains satisfied (assuming = Ufc) after and Qk are changed. The construction of P' and Q' 
is simply to satisfy (|22]i-(|24li. Since 

The construction of P', Q' and can be done recursively. Through its construction given in Algorithm [T| w'{() 
may only violate pT] ) and (|30]l of all the constraints in RMOPF. ^ 

B. The Proof of Main Result 

The first question is whether w'(e) is feasible for RMOPF-line. It turns out that, if the power flow Pm-i, Qm-i, 
. . ., Pq, Qo increases, then not only t(;'(e) is feasible for MOPF-line, but also w'(e) has a smaller objective value 
than w. 

Lemma 2: The new point w'{e) is feasible for RMOPF-line, and has a smaller objective value than w, i.e.. 



provided that 




N 



N 



1=1 



P'k>Pk, Qk>Qk, fc-0,...,m-l. 



J 

Proof: In Appendix [B 

We illustrate Lemma [2] through Fig.js] In constructing the new point w', we do not change P„i, Qm but reduce t 



(31) 



We illustrate L 



Prr 



m, — 1 



m 



N 



Fig. 5. Illustration of j3T}. 



which is proportional to the real and reactive power loss on (m— 1, to). Hence, P^^i > Pm-i and Q'^^i > Qm-i- 
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Algorithm 1 Primal Construction 



Input: 

{p, q, P, Q, i, v), feasible for RMOPF-line but violates ([26] 

i > 1 \ ii > ' ^ ' >; 

Output: 

w'{e)^{p',q',P',Q',i',v'). 

1) Construct p' and q': 

p' ^p, q' ^ q. 

2) Construct P\ Q' and f: 

• for fc = ?7i + 1, . . . , A^: 

K^Pk, Qk^Qk, e'k^ik-, 

• for k = m: 

• for fc = m — 1, m — 2, 



^ P/. 



fc+i 



rk+ii'k+i +Pk^ 



Q'k ^ Q'k+i - xk+ii'u+i + 9fe; 

. for fc = 0: 

Qfc ^ Q'k+i - ^k+i^'k+i + 9fc- 

3) Construct v': 



for k^l,2,...,N. 




4 ^ K-i + ^{r^Pk + ^kQ'u) - {rl + xl)£'^ 



Intuitively, after increasing Pm-i and Qm-i, power flow P^ and Qj., for fc = m — 2. . . . , 1, 0, should also increase. 
Lemma [2] says that if these power flows indeed increase, then w' is feasible and "better" than w (in the sense that 
w' has a smaller objective value than w). It follows that w cannot be optimal for RMOPF-line. Since w is chosen 
as an arbitrary feasible point of RMOPF-line that violates (|26l), RMOPF-line is exact. 
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Inequality ( |3T] l guarantees the exactness of RMOPF-line, and we now investigate when ([31) holds. With the 
simplified notation for one-line networks, equations ([T4ji-([T5|) and ^T9\ take the form 



qk, 



n>k 



for k & Af. We also have 



prcv 
^k 



n 

fe=i 

i-l 



max{P^'",0} , Q 



fc=i 



for i e TV. 

Lemma 3: Inequality holds, if 



"■11 
-a' 




21 



/or i = 1, . . . , TV, and e /s sufficiently small. 

Proof: In Appendix [C] -^^^ 
Condition ( |32] l is nothing but pO] ) in the one-line network. We combine Lemma |2] and [3] to get the following 



theorem. 



.,TV. 



Theorem 2: RMOPF-line is exact if ( [32[ i holds for i = I, 

Proof: If an optimal solution i(;°p* to RMOPF-line violates ( p6] l, then we can construct a new point w(e) from 
w°P' according to Algorithm [l] If ( [32| ) holds for i = 1, . . . , TV, then we can pick sufficiently small e such that 
w{() is feasible for RMOPF-line and has a smaller objective value than w°p', according to Lemma |2]and[3] This 
contradicts with the assumption that u;°p* is optimal. Hence, any optimal solution u'°p* to RMOPF-line must satisfy 
( |26l l, therefore RMOPF-line is exact. ■ 

Generalization of Theorem 12] in tree networks is Theorem fTl whose proof is provided in Appendix iDl 

C. Connection with Prior Work 

In ifTTll . four sufficient conditions for the exactness of a convex relaxation for OPF are given, using the dual 
method. There are mainly two differences between the convex relaxation in ifTTl and RMOPF. Firstly, the relaxation 
in fnl is 



RMOPF is I - ^ ' P2 



> [Pfj + Qfj)/vi with Pij and Qij denoting the power flow from i to j, while the relaxation in 
> {Pj^ + Q'ji)/vj with Pji and Qji denoting the power flow from j to i. Secondly, the voltage 
constraint in ifTTl is w ^ w; while the voltage constraint in MOPF is ^; ^ w, i;'™ ^ v. The four sufficient conditions 
in lEl are 



(pi) Plf^m > 0, Q'-(g) > for e C. 
(p2) r,j/xij = rki/xki for any {k,l) e C. 
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(p3) rjj/x.j < rjk/xjk for any (j, k) e C, and 

i^j"(p)>Ofor(j,j)e£. 
(p4) r.j/xij > rjk/x.jk for any (i, j), (j, fc) e £, and 

gi^"(g)>Ofor(z,j)e/:. ^ 
We establish similar sufficient conditions for the exactness of RMOPF in the following theorem, using the primal 
method. 

Theorem 3: RMOPF is exact if any of the following conditions holds. 
(ci) = 0, Qjr® = 0/or j) e C. 

(c2) r.,.j/xij = rfcz/xfc; /or any (i, j), (fc, /) € £, and 

~ 2r,,i-r(p) - 2x„Q5r(g) > 0/or j) G C. 
(c3) r^j/xij < rjk/xjkfor any {j,k) G £, and 

(c4) r^j/xij > rjk/xjk for any {j, k) G C, and 

Q]Ti<l) = 0, - 2r„PjT(p) > 0/or G £. 
Proof In Appendix |E] 

We compare conditions pl-p4 with conditions cl-c4. Firstly, they can all be checked before solving the corre- 
sponding relaxation. Secondly, since P}™ — ~Pj"^ and Q\"^ ~ ^Q^ji ^'^^ ^ conditions pi and cl are 
identical. Thirdly, Vj — 2rijPjf^ (p) — 2xijQ'jl^ (q) > always holds in practical networks, for every G C. 
Hence, conditions p2 and c2 are identical for practical networks. Similarly, conditions p3 and c3 are identical for 
practical networks, and conditions p4 and c4 are identical for practical networks. 

To summarize, we have established sufficient conditions for the exactness of RMOPF using the primal method. 
These conditions are are identical to the conditions given in ifTTl in practical networks. 





D. Uniqueness of MOPF Solutions ^ 

We study the uniqueness of MOPF solutions in this section. We assume that RMOPF is exact, which implies a 
point solves MOPF if and only if it solves RMOPF. Therefore we only need to study the uniqueness of RMOPF 
solutions. When the objective function of RMOPF is strictly convex, the solution is unique. However, the objective 
function is often linear — not strictly convex, e.g., the power loss. We show that with arbitrary convex objective 
functions (not necessarily strictly convex), RMOPF has a unique solution if it is exact. 

Theorem 4: If RMOPF is exact and its objective function is convex, then RMOPF has a unique solution. 
Furthermore, MOPF has the same unique solution. 

Proof: In Appendix |F] ■ 

V. Conclusion 

We have proposed a convex relaxation for solving a modified optimal power flow problem in tree networks. A 
sufficient condition for the exactness of the relaxation is given. This condition can be checked before solving the 
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relaxation, and holds in various test networks. 
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Appendix 

A. Proof of Lemma [7] 

Since {p,q, P,Q,£,v) satisfies Q-Q and we have i.j > [Pj,^ + Q]^) jvj > for all € C. It follows 



k-.jGVk 

lin 



for e C. Similarly, Q,, < Qf^{q) for e C. 
It follows from (|4l) that 



' V 




Vj Vi — 2{rijPji + XijQji) i^ij ^ 
• < 2(r,,Pjr+x,,C^)^ 
for e C. We sum up the inequalities over Vj to obtain 

vj-vo< E 2irkiPl'l^' + XkiQfn^vT\P,q)^vo, 

which implies Vj < q), for j = 1, . . . , iV. 

Vr 

By the construction in Algorithm [T] w' satisfies (|22]l-(p5]l and (|28|)-(|29ll. To show that w' is feasible for RMOPF- 
hne, we are left to check that w' satisfies ( p7] i and ([30|. 
Claim 1: The voltage v' y v if ( |3T] ) holds. 

Proof Define AP := P' - P, AQ := Q' - Q and Av := v' - -y. If ([31) holds, then APfe > 0, AQk > for 
= 0, . . . , m - 1, and APk = 0, AQk = for fc = m, . . . , iV. It follows from (|23|-(|25| that 



Vi = v^-i + ri{Pi + Pi^i -pi-i) + Xi{Q., + Qi^i - qi-i), 

which implies 

Avi = Avi-i + n{APi + AP.^i) + x^{AQi + AQi^i) 
> Av.^u 
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for i — 1, . . . , N. Moreover, the inequality becomes strict when i = I since APq > 0, AQq > 0. Hence, 

Aujv > • • • > Av2 > Avi > Avo = 0, 

which completes the proof of Claim [T] 

It follows from Claim [T] that v'^ > Vi > and 

^, ^ max{J^'^J^n + niax{Q,f,Qf } 

> Pl + Ql^Pl + Ql 

Vi v[ 

for i = 1, . . . , iV. Hence, w' satisfies ( p7] i and ([30|, therefore feasible for RMOPF-line. It has a smaller objective 
value than w because 

N N 




Y.-A-Y.rA ^^^^ 



i=l i=l 



Vi=0 1=1 / \i=0 i=l / 



Pn - ^0 > 0. 




C. Proof of Lemma |5] 

Lemma 4: Given (ai, . . . , a„), . . . , 6„), (ci, . . . , c„), (di, . . . , (i„) satisfying < a/c < 1, 6^ > 0, > 0, 
< rffc < 1 /or fc = 1, . . . , n, where n € {1, 2, . . .}, define matrices / 

ykA. = 1 I (33) 



-Ck 1 - dk 



for k = 1, . . . , n. Let Un+i y be a 2 x I vector, and define 
for k = l,...,n. If 




71 n 

y^cfc j]^ (1 ~ 4 



Un+l >- 0, 



k=l fc=l 




f/ien Ufc ;^ /or fc = 1, . . . , n + 1. 

Proof: We prove Lemma |4] by induction on 
(i) When n = 1, Lemma [4] is trivial. 
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(ii) Suppose Lemma |4] holds for n = 1, . . . ,K (K > 1). When n = + 1, we have 

/ n n \ 

k=2 k=2 



V 



k=2 



T{{l-dk) 



Un+l 



k=2 



/ 71 n \ 

k=l k = l 



Un+l >- 0. 




\ fe=i fe=i / 

According to induction hypothesis for n — K , Uk for k — 2, . . . ,n + 1. To complete the induction, we 
are left to show that ui >- 0. The idea is to construct a new sequence of K matrices and apply the induction 
hypothesis foin — K again. The construction is 



1- (l-ai)(l-a2), b[ := 61+62, 
ci + C2, d[ := 1 - (1 - di){l - (£2); 



afc+i, 6'^ := 6fc+i, := Ck+i, and d'j, d^+i for k ^ 2, . . . , K. Define matrices 



A' ■— 

^k • — 



-c', 1 - d[. 



for fc = 1, . . . , and uj, := A'^ ■ ■ ■ A'j^Un+i for fc = 1, . . . , iiT. Apply the induction hypothesis for n = K, 
and we obtain u'^ y 0, u'2 y 0. Then 

= AiA2A3...AK+lUn+i 
= AiA2^2 • ■ ■ A'l^Un+i 
= AiA2u'2 h A[u'2 
= m'i >- 0, 

which completes the induction. 
According to (i) and (ii). Lemma |4] holds for n ~ 1,2, . . .. ■ 
^ Recall that we are given a point w = [p, q, P, Q, £, v) that is feasible for RMOPF-Une but violates ( |26] l, and m 




is the first bus that violates (l26|. Given any e satisfying < e < 



£,n, we can construct a point w'(e). 
What we want to prove in Lemma [3] is that if inequality ( [32] i holds and e is sufficiently small, then inequality pT| ) 
holds.^ 

Define non-negative numbers 



2rfePr(p) 


j: 2r,Or(9) 
Ok ■ — 


Vk 


V-k 




2xkQr{q) 

«fc •= 


V-k 
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for k = 1, . . . , N . Inequality p2| ) is equivalent to 

a\iri > a\2Xi > and 022^;^ > a2iri > 0, 

which implies a\i > and 022 > 0, for i = 1, . . . ,N. Since a\i and 0^2 are products of 1 — aj; and 1 — dk 
respectively, we have Ofc < 1 and dk < 1 for fc = 1, . . . , iV — 1. For a; e M, define ma.x{x, 0}. Define 



Ck := 





_ 2rkQt 


, Ok ■ 




Vk 


Vk 


2xkP^ 

, elk 




Vk 


Vk 



for k = 1, . . . , N, then we have 



< Ofc < afe < 1, 0<bk<bk, 



for fc = 1, 



, AT - 1. It follows from ([32li that 



>- 



0<cfc 


< Cfc, 


< dfe < dk 


32| that 






m — 1 

m— 1 


Ofc) 


m-1 \ 

m—1 




Cfc 


n (i-^'') 
fc=i / 


m— 1 


Sfc) 


m-1 \ 

m—1 


fc=l 




n (1-^0 

i:=l / 



2> 



X7J 



Xr. 



'ii 



*12 



m m 
"21 "22 

Define matrices Ak according to ( [33] l for = 1 , . . . , A^, and 

^m ■ — (^m;*^m) 

• • j4,„_iu,„ >- for fc = 1 



1 









^ 0. 



It follows from Lemma |4] that Ak 




^Pk 
AQk 



, m. We will show that this implies 

= Ak+i ■ ■ ■ A„^-iu,ne + 0{e^) y 



(34) 



for = 0, . . . , m — 1, if e is sufficiently small. 
We prove ( |34) l by induction on k. 
i) When k = m — 1, it follows from Algorithm 1 that 

APk 
^Qk 



e = w,„e >- 0. 
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ii) Suppose ( [34| l holds for k = t (1 < t < m — 1), we will show that ( (34[ l also holds for k = t — 1. We have 

max{P/2 , P2 } + niax{Qf , } + q2 



nmx{P/2 - p2, 0} + max{Qf - Q?, 0} 
max{2PtAPt + O(e2),0} 



+ 



max{2QtAQt +O(e2),0} 



2P+APt , 2Q+AQt 



when e is sufficiently small such that APf > and AQt > 0. It follows that 

APt_i \ / APi 

Agt_i ) \ AQ, 

V AQ 

= At • • -^m-lWrne + O(e^). 

Since Af ■ ■ Am-iu,n >- 0, ([34]) holds for fc = t — 1 if e is sufficiently small. 
According to (i) and (ii), ( (34] i holds for fc = 0,...,TO— lifeis sufficiently small, which completes the proof of 
Lemma [3] 




D. Proof of Theorem |7] 

If a solution w°p' = {p°P\q°P\ P°p\Q°p\ £°p\v°p^) to RMOPF violates ^, then there exists a line (i, j) such 
that 




tj opt 

and (jS]) holds on Vi- For e satisfying ^ 

J w < . < i^rijiMn! ,opt 

we construct a new point 'w{e) by 

1) setting p = p°P*, g = q°P^; 

2) setting 4; - C'' A.; = P^r- - QT (A:, ^ P,; 

3) setting f^;, Pki, Qu as Algorithm 1 for {k,l) g P^; 

4) setting v to satisfy Q. 
Define 

Po := Po + X! ~ ''ofc^fc) , 

(o,fc)e£ 

Qo '■— QO + ^ (QfeO — Xok^ok) 

(o,fc)e£ 
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as the power flow from the substation to the main grid. Follow Appendix |Cj we can prove that if ( |20] l holds, then 
the point w(e) satisfies 

APik > 0, AQik > 

for {k, I) £ Vi and APp > 0, AQq > for sufficiently small e. Follow Appendix [B] we can further prove that such 
w(e) is feasible for RMOPF and has a smaller objective value than w°p*, which contradicts with the assumption 
that ti;°P* is optimal for RMOPF. Hence, any optimal solution w°p' to RMOPF satisfies Q, therefore RMOPF is 
exact. 




E. Proof of Theorem [i] 

We only give the proof for a one-line network as in Fig. |4] The proof can be generalized to tree networks 
following Appendix |D] With the simplified notation in a one-line network, we re-state Theorem [3] as 
Theorem 5: RMOPF-line is exact if any of the following conditions holds. 
. ^p) = 0, = /or A: = 1, . . . , TV. 

• r„i/xm — Tn/xn for any 1 < m,n < N, and 

- 2rkPrip) - ^XkQriq) > 0/or fc = 1, . , 

• Tm/xm < fn/xn for 1 < m < 71 < N, and 

prm = 0, - 2xkQr{q) > o/or ^ = i, . . 

• Tm/xm > Tnjxn for 1 < m < 11 < N, and 
Qriq) =0,Vk~ 2rkPrip) >Ofork=l,.. 

We prove Theorem |5] through Claim |2||5] 

Claim 2: RMOPF-line is exact i/Pf ^(p) = 0, Q't^^iq) ^ for k = 1, . . . , N. 

Proof If Pf ^(p) = and Ql?''iq) = for k = 1, . . . , N, then a\i = a^^ = 1, al^ = for i = 1, . . . , N. 

It follows that ([32]) is satisfied for i = 1, . . . , iV, therefore RMOPF-line is exact according to Theorem [2] ■ 

Claim 3: RMOPF-line is exact if Tm/xm — Tn/xn for 1 < m,n < N, and Vf. — 2rkP^^{ji) — 2xkQ]?''{q) > 
fork = l,...,N. 

Proof: Let 77 := ri/xi denote the r/x ratio for all transmission lines. We will show that ( (3T| holds for 
sufficiently small e. Since 



JV 




AP,. = - ^nM„ 



n=k + l 
N 



= -77 ^ a;„A^„ = -riAQk 

n=k+l 

for fc = 0, . . . , iV, it suffices to show AQ^ > for fc = 0, . . . , m — 1 to prove ( |3T] i. For brevity, we ignore the second 
order term O(e^) in the following derivation, but the second order term can be considered following Appendix [c] 
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It follows from (|35) that 




= - (wt - 2nP+ - 2xtQt) ^Qt 

Vt 

for i — 1, . . . , m — 1. Therefore AQf has the same sign for t — 0, . . . , m — 1 since vt — 2rtP^ — '2xtQf >1 
- 2rtPl°"'{p) - 2xtQ\°"'{q) > for t = 1, . . . , m - 1. By construction, A(3,„_i = a;„ie > 0. Hence, AQt > 

for t = 0, . . . , m — 1. Therefore ( (3T] i holds, which implies exact RMOPF-line. iT'^^^^Ilk ^"'^ ' 

Claim 4: RMOPF-line is exact ifr,n/xrn < rn/xnfor \<m<n<N,and Pl°^{p) = 0, v^.-2xkQ]^'' {q) > 

Proof: We will show that ( (3T| holds for sufficiently small e. For brevity, we ignore the second order term 
O(e^) in the following derivation, but the second order term can be considered following Appendix [c] 
We have < P^t < P'k^{p) = 0, which impUes P+ = 0, for A: = 1, . . . , A^. It follows from ([35) that 

AP,_i \ ^ \ I 

AQt^i ) \0 I- ) \ AQt 

for t — 1, . . . , m— 1. Therefore AQt has the same sign for t = 0, . . . , m— 1 since l~2xtQf /vt = ^ {vt ~ 2xt(5"^(q)) > 
for t 1, . . . , m — 1. By construction, AQm-i = > 0. Hence, AQt > for t — 0, . . . ,m — 1. Now we 
prove that _ 



An ^ (36) 

AQt Xm 



for i = 0, . . . , m — 1 by induction on t. 
i) When t — rn — 1, we have 



^ APt ^ r.,^ W 



APt _ rvj£ _ Tn 

InequaUty ([36| holds, 
ii) Suppose (l36| holds for f = /s (1 < fc < m — 1), then 




AP,_i AP, - ^AQ, 



AQfc_i (i-^ 



1 + 



> 



2 _ a;,^ 



i.e., ( [36| ) also holds for i = /c — 1. Here the first inequality is because ( [36| l holds for t = fc (induction hypothesis); 
and the second inequality is because r„xlx„x > Vk/xk- 
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According to (i) and (ii), ( (36[ l holds for < = 0, . . . , m — 1. It follows that APt > for < = 0, . . . , m — 1. We have 
proven that ( (3T] i holds, which implies exact RMOPF-line. ■ 

Claim 5: RMOPF-line is exact ifr^^jx^n < r„/a;„/or 1 < m < n < N, and Q]?^{q) = 0, Vk-^rkP^^ {q) > 
/or fc = 1,...,7V. 

Proof: The proof of Claim |5] is similar to that of Claim |4] and omitted for brevity, 

F. Proof of Theorem [4] (^^j^ 
Claim 6: RMOPF has a unique solution if it is exact and its objective function is convex. 

Proof Let / denote the convex objective function of RMOPF, let w — (ji, q, P, Q,i,v^, w = (ji, q, P, Q, £, 
denote two arbitrary solutions to RMOPF, and defnie 

as the optimal value of RMOPF. Since RMOPF is exact, 

^ij^j ~ -^ji ^" Q ji^ ij^j ^ ji ~^ Q ji ^ '^^^ (37) 

for E C. For any 9 E (0, 1), define 

w{e) -.^Ow + ii-e)!!) 

The point w{d) is feasible for RMOPF since RMOPF is convex. It is also optimal for RMOPF because 

f{w{e))<ef {w) + {i-e)f{w) = f°'p\ 

If RMOPF is exact, it follows that ^ '"^^ 

e,,{9)v,{e) = pf,{e) + Q%{e) 08) 

for e C. Substituted^ * ^^^Ib^ 

p i,j{0) = ^4 + (1 - 6%,, v,{9) = + (1 - 9)vj, 
A P^{9) = 9P,j + (1 - 0)4, Qj,{9) = 9Qj, + (1 - 9)Qj, 
into ( |38] l and simplify using ( [37| ) to obtain 

for e C. It follows that 

"''^^ 2 ^Pj^Pj^ -\- QjiQji 

p2 , /)2 p2 , q2 

P2. p2 02. 02 

= + V. + V. + W," 

Vj Vj Vj Vj 




> 2^Pf^P^ + 2^Q%Q 



^ ^Pji^ji ^" 2QjiQji. 
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for E C. The two inequalities must both attain equalities. The first inequality attaining equaUty implies that 



p 




Qji 






Qji 



for £ C. The second inequaUty attaining equality implies that Pji and Pji have the same sign, Qji and Qji 
have the same sign, for G C. It follows that 



P,: 



Q, 



for G £, which further implies 



for e C. Define 



for j £ Af, then 770 = 1, and 



P. 



P. 



31 



Q^ 

Qji 



o 



Vj 2{rijPji + XijQji) + (rjj + Xij)iij 



= £, and v = v. It is not difficult 



Vj li^TijPji -\- XijQji^ -\- {j^ij ^" "^ij^^ij 

for {i,j) G C It follows that rjj = 1 for j E Af, which implies P = P, Q = Q 
to further prove that p = p and q — q, which completes the proof of w = w. Claim |6] follows from w = w. ■ 
To complete the proof of Theorem]?] we only need to show that RMOPF and MOPF has the same solution when 
RMOPF is exact. Any solution w' to RMOPF is also optimal for MOPF because RMOPF is exact. For any solution 
w°P' to MOPF, we have f{w°P^) < /(w'). Hence, w"^^ is optimal for RMOPF. We have proved that solutions to 
MOPF and RMOPF are the same when RMOPF is exact, and Theorem |4] follows. 
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